Damping and coherence in a high-density magnon gas 
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We report on the fast relaxation behavior of a high-density magnon gas created by a parametric amplification 
process. The magnon gas is probed using the technique of spin-wave packet recovery by parallel parametric 
pumping. Experimental results show a damping behavior which is in disagreement with both the standard model 
of exponential decay and with earlier observations of non-linear damping. In particular, the inherent magnon 
damping is found to depend upon the presence of the parametric pumping field. A phenomenological model 
which accounts for the dephasing of the earlier injected magnons is in good agreement with the experimental 
data. 
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I. Introduction 

Several decades after the first phenomenological descrip- 
tions, the nature of spin-wave dissipation is still a subject of 
intense scientific research 1-6 . An understanding of spin- wave 
damping mechanisms is crucial both for the furtherance of the 
field of fundamental spin dynamics 7 and practical applications 
in novel magnetic memory, microwave, and spin-wave logic 
devices 8 . Recent developments in the study of energy transfer 
in magnon systems, including the discovery of Bose-Einstein 
condensation of magnons 9 , have attracted particular interest to 
the problem of nonlinear damping and decoherence of densely 
populated magnon states, created by the effect of parametric 
pumping 1 2il!. 

Up to now there has been no access to the information 
on relaxation of highly populated short- wavelength magnon 
states as these spin waves are too short to be directly de- 
tected by conventional microwave and optical methods. Only 
measurements of thresholds of parametric instability provide 
us with information about the decay characteristics of these 
magnons at low population level. However, even in this case 
the magnons are influenced by an external microwave pump 
field. To overcome this obstacle and to access the freely evolv- 
ing highly populated magnon groups, we used the recently de- 
veloped technique of spin- wave signal recovery 12 . 

In this work we demonstrate that free evolution of a dense 
group of parametrically injected magnons can be described 
neither by its monotonic exponential relaxation nor by possi- 
ble nonlinear decay. The phase decoherence of the frequency 
smeared magnon group should be taken into account in or- 
der to understand the experimental data. We chose a single 
crystal yttrium-iron-garnet (YIG) ferrimagnetic film as a test 
object in our experiments. Due to its extremely low natural 
magnetic damping 13 this material is widely used in works on 
magnon gases and condensates 7 . Furthermore, currently one 
observes a clear revival of interest to YIG films due to the re- 
cent discovery of the spin Hall 14 and spin Seebeck 15 effects 
in Platinum/YIG bi-layers. Moreover, submicron-thick poly- 
crystalline YIG films are now promising as candidates for 
novel microwave devices 16 , and conventional single-crystal 
micron-thick films have been recently used to demonstate mi- 
crowave encoding 17 . 
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FIG. 1 . (Color online) a) Experimental setup, consisting of a single- 
crystal YIG waveguide, input and output antennas and a dielectric 
resonator for the application of the pumping microwave field. The 
waveforms are the shapes of the input, pumping and restored pulses 
as seen on the oscilloscope, b) Section of the spin- wave dispersion 
spectrum for a thin magnetic waveguide schematically showing the 
hybridization of Damon-Eshbach (dashed line) wave with standing 
spin wave modes. 



II Experiment 

A 100 ns-long traveling spin- wave packet (ri npu t = 100 ns) 
is excited at a frequency of 7 GHz by a microstrip input an- 
tenna in a 1.5 mm- wide and 5 um-thick YIG stripe. An exter- 
nal magnetic field H of 143.64 kA/m (1805 Oe) is applied in 
the plane of the stripe, perpendicular to the spin- wave propa- 
gation (Fig.Hk), i.e. in the Damon-Eshbach (DE) geometry 18 . 
The excited spin- wave packet traverses the magnetic sample, 
is detected by the output antenna placed 12 mm apart from the 
input one, amplified, and observed with an oscilloscope. After 
450 ns since the passage by the wave packet of the central area 
of the YIG stripe (to = 450 ns), a pumping microwave field 
h p at 14 GHz produced by a dielectric resonator is applied par- 
allel to the bias magnetic field. As a result of the pump action, 
an additional "recovered" pulse is observed at the output an- 
tenna at a frequency of 7 GHz (see output waveform in Fig.Hk 
and dashed line in Fig. [2]). Since the DE packet has already left 
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FIG. 2. (Color online) Spin- wave power envelope as experimentally 
observed at the oscilloscope. The light red areas schematically show 
the time intervals when the pump is on. The continuous red line 
represents the results of the two pump-pulse recovery. The first pump 
pulse is switched off before the recovered signal saturates. After 
a delay time Tdeiay the pumping is switched on a second time the 
second recovered signal is observed. The dashed blue line represents 
the experimental result for the pumping not being interrupted. 



the area of parametric interaction, this is not the result of the 
direct amplification of the signal puls o 19 i 20 . Rather, the para- 
metric pumping acts on the standing spin- wave (SSW) modes 
(see Fig. [TJ)) which are excited across the film thickness by 
the traveling wave packet through a two-magnon scattering 
mechanism. These modes are then amplified by parametric 
pumping and are scattered back to form a new traveling DE 
wave, which is finally detected at an output antenna in a form 
of the "recovered" pulse (Fig.[T^). 

Unlike in previous studie s 12 ! 21 ! 22 , in the experiment reported 
here, two consecutive parametric pump pulses are used rather 
than one. The response of the magnonic system to the second 
pulse is utilized as a probing tool for extracting information 
about the free relaxation behavior of parametrically excited 
magnons during the pump free pause. Therefore, we cut off 
the first recovered pulse by stopping the parametric pumping. 
After a certain delay time Tdeiay, a second pumping pulse is 
applied and a second recovered pulse is observed (see contin- 
uous line in Fig. [2]). We register the peak amplitude A max and 
the arrival time t max of the second recovered pulse. 

As one sees from Fig. [3j the behavior of the second re- 
covered pulse strongly depends on the duration of the first 
pumping pulse T pump i for small T pump i < 450 ns as well 
as on the delay time Tdeiay between two consecutive pump- 
ing pulses. It is remarkable that the latter dependence is non- 
monotonic: with an increase in Tdeiay the second recovered 
pulse first grows, reaches a maximum, and then decreases in 
intensity (Fig. [3^). In particular, one observes a gain in the 
amplitude of the second recovered pulse compared to the ex- 
periment in which the pumping is uninterrupted (Tdeiay — ^ 
in Fig. [3k). 

Ill Theoretical model 

One easily finds that neither the initial increase in the am- 
plitude nor the subsequent non-exponential decrease (note the 
logarithmic scale) in the peak amplitude can be explained in 
the framework of the standard approach based on association 
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FIG. 3. (Color online) a) The amplitudes of the second pulse shown 
in Fig. [2k are extracted and plotted against Tdeiay for different du- 
rations Tpumpi of the first pumping pulse. The amplitudes are nor- 
malized to the maximum observable amplitude using only one long 
pumping pulse, b) The time t ma x at which the peak amplitude of the 
second pulse in Fig.[2k is observed as a function of Tdeiay 



of a specific relaxation time with each decay process. For 
this reason we make an attempt to employ a more detailed 
model. The model we use is an extension of the theory in 
Ref. 22 . The latter succeeded in describing the single-pulse 
pumping. Treating a two-pulse pumping with a significant 
length of the pause between the two pulses requires consider- 
ing dephasing within the parametrically excited wave packet 
during the pause. Thus, in contrast to Ref.—, one cannot use 
ensemble-averaged equations for magnon densities and has 
to consider each magnon group as consisting of a number of 
pairs of waves with different eigenfrequencies. 

We start the description of the developed theory by noting 
that the recovered pulse has a pulse-like shape. This sug- 
gests that competition of two magnon groups for energy pro- 
vided by the parametric pumping process 23 takes place in the 
magnon system, since a pulse-like shape of the observed sig- 
nal is a typical result of this competition. This process can 
be explained in the following way. The parametric amplifica- 
tion is frequency selective: only magnon groups with eigen- 
frequencies ujk which lie inside the narrow frequency band 
cjp/2 — v < uo(k) < cjp/2 + v whose width is equal to the 
parametric amplification gain v — h p V are amplified (V is the 
parametric coupling coefficient). They are driven at half the 
frequency u; p /2 of the applied microwave field. Importantly, 
there is no restriction on the magnitude and the direction of the 
magnon wave vector: the process allows creation of magnons 
with arbitrary wave vectors, provided their frequencies are in- 
side this frequency band. Thus, given the large number of 
frequency-degenerate magnon dispersion branches, magnon 
groups over a large range of wave vectors are amplified. This 
includes externally excited oscillations, hereafter called the 
signal group, as well as thermally activated magnons with 
wave vectors up to 10 5 cm -1 . The latter magnons, which 
are considerably decoupled from structural defects of the YIG 
film because of their short wavelength, experience the lowest 
two-magnon decay and consequently show the highest para- 



metric gain 24 . Hence, we refer to this dominantly amplified 
magnon group as the dominant group. 

As the phase shift between the external pumping field h p 
and a microwave magnetic field induced by the parametri- 
cally pumped magnons ("internal pumping") increases with 
increase in the total magnon density, the resulting effective 
pumping field inside the sample consecutively diminishes 22 . 
At some point in time the total magnon density reaches a 
critical threshold level A cr at which the effective pumping 
becomes so small that it is able to further support just one 
magnon group, the one which has the lowest decay rate and 
the highest efficiency of parametric interaction: the domi- 
nant group. This leads to the suppression of the signal group, 
whilst the dominant group reaches saturation at A cr . From this 
time on the system stays in a quasi-equilibrium state wherein 
the magnon density of the dominant group is saturated. The 
signal group is exposed to an effective pumping below the 
generation threshold 10 and thus should decay with time to the 
thermal level. As a result of this process, the signal group, and 
consequently the measured signal, acquire a pulse-like shape 

To obtain a quantitative description of the magnetic dy- 
namics for the present case of two pump pulses, we per- 
form numerical simulations based on the nonlinear model of 
the magnon group competition. We assume that the signal 
group represents a set of spin-wave oscillations (standing- 
wave modes) C^ with resonance frequencies uj uniformly dis- 
tributed across the frequency band Suj = ±1.5 MHz around 
the half-pump frequency uj p /2. This frequency band corre- 
sponds to the spectral width of the dipole-exchange gaps in 
the spin- wave dispersion (Fig. lb). Since the recovered pulse 
is only observed in these gaps, only magnons within this fre- 
quency band should contribute to the recovery process. 

The standing-wave modes are excited through 2-magnon 
scattering processes by the traveling-wave pulse on its way 
across the film. For simplicity we assume that the initial phase 
of magnetization precession for all the modes is the same. The 
output traveling-wave pulse is formed due to 2-magnon back- 
scattering and thus represents an exact time replica of the dy- 
namics of the signal group. Therefore it is sufficient to con- 
sider only the evolution of two immobile magnon groups: the 
signal and the dominant ones (The dominant group consists of 
magnons with complex amplitudes C d .) 

The general dynamic equations are as follows: 

[d/dt ± r s ± i(Co s - u p /2)] ci - iP*c s j = o, (i) 

[d/dt ± T d ± i(cD d - u p /2)] C d - iP^Ct* = 0, (2) 

where the stars denote complex conjugation, Y s and Yd are the 
relaxation parameters for the signal and the dominant group 
respectively, 



uj s = w + Aw ss + Aw sd , 


(3) 


w d = u + Aw dd + Aw ds 


(4) 



are the eigenfrequencies of standing-wave oscillations with an 
account of nonlinear frequency shift for both signal and dom- 
inant groups respectively, and 

P sd , 
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(5) 
(6) 



is the effective pumping for the signal and the dominant 
groups respectively. In Eqs. (3-6) and below the upper in- 
dices "s" and "d" denote the signal and the dominant group 
respectively. Double indices of type a/3 denote the action of 
the group /3 on the group a. If both indices are the same the 
respective coefficient or magnitude describes self-action. 

In our calculations we make a number of important simpli- 
fications which allow us to considerably decrease the compu- 
tation time. First we use a usual assumption 22 that the non- 
linear frequency shift is not important for the dynamics of the 
dominant group (Aw dd = 0). Second, we do not need to take 
into account the contribution to the total nonlinear frequency 
shift by the signal group (Auj ss = Auj ds = 0), since its am- 
plitude is considerably smaller than the amplitude of the dom- 
inant group at the time, when the effective pumping saturates. 
Thus, the only term we have to account for is 



Auj sd 



T sd J2 C d ,C d t,C^,5(uj'-uj"+uj'"-uj), (7) 



where S denotes the Kronecker Delta, and T sd is the respec- 
tive nonlinear coefficient. In the following we will use a short 
hand notation T sd = T. 

Similar considerations apply to the contributions to the to- 
tal pumping which consists of the external pumping (y) and 
the internal pumping given by the remainder of the terms in 
Eqs. (5, 6). We obtain: P ss = P ds = 0, 



psd = ^sd J2 C d ,C d „C s J„, S(uj f 
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ydd V^ 
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Ct>Ct»C%» S(w' - w" + J" - w),(9) 



where S afS are the respective nonlinear ff-coefficient a 22 ! 23 . 

Test numerical calculations show that calculation results do 
not qualitatively change with the variation in the difference 
S dd — S sd in reasonable limits. Therefore, to minimize the 
number of degrees of freedom we set S dd = S sd = S and, as 
previously, we assume that S > 22 . 

The last simplification we make is removing the frequency- 
mixing terms from Eqs. (7-9). This simplification does not 
lead to a qualitative change in the results of our numerical cal- 
culations, but enormously decreases the computation time. It 
reduces Eq. (4) to a simple formula for the effective pumping: 



P s , = p: 



d _ 



Pu, = v + S^2(C*,) 



(10) 



For the pump-free period the same equations (1-10) are 
valid, one just assumes v = 0, so that only the internal pump- 
ing is on during the pause. 

The structure of obtained equations allows us to renormal- 
ize the spin-wave amplitudes G u such that S = 1 and only 
the ratio \S/T\ matters (T < 0). Initial conditions for the 
dominant group are random distribution of amplitudes for its 
frequency components. Initial conditions for the signal group 
follow from the expression for the frequency spectrum of a 
rectangular pulse: 



C w {t\) = Coex.p(iujti)F(Suj) 
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FIG. 4. (Color online) a) Left axis: measured (circles connected by 
a dashed line) and simulated (dashed line) time t max at which the 
peak amplitude of the recovered pulse is observed as a function of 
Tdeiay for T pU mpi = 400 ns. Right axis: measured (dots connected 
by a solid line) and simulated (solid line) peak amplitudes of the 
recovered pulse for T pump i = 400 ns. b-d) Simulated time profiles of 
the signal (thin solid lines) and of the internal pumping (bold dashed 
lines). The light red areas indicate the time intervals when the pump 
is on. The gray area in panel b) indicates the input spin- wave pulse, 
b) Uninterrupted pumping (rdelay = 0). c) Tdeiay = 300 ns. d) 
Tdeiay = 700 ns. t = corresponds to switching on the second 
pump. Duration of the first pump pulse T pU mpi = 400 ns. 



where T tw is the length of the initial traveling wave pulse 
which generates the standing wave oscillations through the 
2-magnon scattering process, and To is the time interval be- 
tween the the passage of the traveling wave pulse through 
the pump area and the begin of the first pump pulse, Co is 
some constant, and F(Suj) is equal to one within the dipole 
gap in the dipole exchange spectrum of travelling spin waves 



(Fig. 1), where the parametric signal recovery is usually ob- 
served (F(Suj) = 1, for —Suj/2 < uj < Suj/2) and vanishes 
outside the gap (F(Suj) = 0, for \lj\ > Suj/2). 

The system of equations (1,2) was solved numerically. We 
used 200 discrete values of eigenfrequency detuning uj from 
the half-pump frequency in the range from —2v to 2v. This 
formed a system of 200 coupled nonlinear equations which 
was resolved by using 4th-order Runge-Kutta method. 

We do numerical calculations for r pump i = 400 ns and for 
a number of ratios | S/T | , a number of values of rates of para- 
metric gain, and a number of decay rates for both magnon 
groups. Each time all these parameters, including Co, are 
chosen such that the restored pulse peaks at the same time 
as in the experiment when the parametric pumping is not in- 
terrupted (rdelay = 0). A natural constraint we impose on the 
validity of simulation results is a single-peaked shape for the 
restored pulse for all experimental lengths of the pump-free 

pause Tdeiay 

IV Discussion 

A number of simulation runs allows us to locate the area in 
the parameter space, where the simulated behavior is close to 
the experimental one displayed in Fig. 3. Figure 4a shows the 
best fit we obtain. One sees a fair agreement with the exper- 
iment of the simulated peak amplitudes of the restored pulse 
and of the times for the peak arrival t max as a function of the 
pause length Tdeiay This calculation shows that the terms in- 
volving S and T coefficients do not contribute to the dynamics 
during the pump-free pause, as the parametric amplification is 
quite far from saturation during this time interval. Note that 
analytical solutions exist for S = T = and they are in full 
agreement with our simulations for the first two stages of the 
considered process (t < in Fig. 4b-4d). 

Figures 4b-4d show the behavior of the macroscopic ampli- 
tude of the signal group A s = 10 log(| J^ ^| 2 ) an d of the 
strength of the internal pumping A^ = 10 log(| ^2 U , C^, | 2 ). 
One sees that during the pump pause both magnitudes do not 
vary linearly on the logarithmic scale. The signal-group be- 
havior is close to parabolic on this scale. Note the increase in 
the amplitude for the signal group during the first 200 ns after 
switching off the first pump in Fig. 4d. A careful analysis of 
this stage shows that the increase is due to phase reversal in 
a parametric-echo-like process 25,26 . Recall that the coherence 
of the originally deterministic signal is lost due to dephasing 
during the time interval To between the passage of the travel- 
ing wave pulse and the application of the first pump pulse. In 
the phase reversal process the coherence is restored. The full 
phase restoration occurs for Tdeiay = to- 

The behavior of Ad is more complicated: for the first 150 ns 
after the first pumping has been switched off, the internal 
pump decreases parabolically on the logarithmic scale, then 
its behavior switches to more or less linear on the same scale. 
This non-monotonic behavior of the internal pump also orig- 
inates from dephasing. During the first pump pulse the fre- 
quency width of the dominant group reduces significantly due 
to the frequency selective amplification of the initially thermal 
signal. The phases of the waves which belong to this group 
are locked to the external pumping and the behavior of the in- 



ternal pumping follows the linear behavior of the number of 
magnons, which is in full agreement with L'vov's S-theory. 
However, once the first pumping has been switched off, the 
phase coherence within the group is lost and the magnitude of 
Ad starts to decrease. As a result, at large times (t > 300 ns) 
the magnitude of Ad essentially follows the exponential de- 
crease in the number of magnons, but deviates significantly 
from this law because of the dephasing for t < 300 ns. 

Both echo-like behavior of the signal group and dephasing 
in the dominant group contribute to the non-monotonic behav- 
ior of the peak amplitude of the restored signal. From Fig. 4 
one clearly sees that for Tdeiay = 300 ns (Fig. 4c) the am- 
plitude of the signal at the end of the pause (t = in this 
graph) is quite close to the amplitude of this group at t = 
for the uninterrupted pump. However, the internal pump by 
the dominant group at t = for Tdeiay = 300 ns is smaller 
than for Tdeiay = by 25 dB. This allows the recovered signal 
developing a larger peak amplitude for Tdeiay = 300 ns than 
for Tdeiay = 0. Obviously, a larger Tdeiay requires a larger 
time for the restored signal to peak, which results in a linear 
dependence of t max on Tdeiay in Fig. 3a. 

In agreement with experiment, the model shows a vari- 
ation of A max and t max with r pump i. However, in the 
model the effect is not so pronounced as in the experiment. 
We also checked if nonlinear damping in the form of four- 
magnon scattering processes involving all spin waves exist- 
ing at the half-pump frequency can contribute to the non- 
monotonic response. Our calculations for T pump i = 400 ns 
have shown that the nonlinear damping is negligible during 
the pump pause for reasonable values of the nonlinear damp- 
ing coefficient 27 . The amplitudes of both signal and dominant 
groups are just too small during the pause to develop the non- 
linear damping. Accounting for the nonlinear damping during 
the second pump pulse does not lead to qualitative changes in 
the behavior either, unless one assumes that the coefficient of 
nonlinear damping is larger than T by several orders of mag- 
nitude, which is unreasonable. 

As a final note of this section we want to comment on appli- 



cability of the two-pulse method for other materials. YIG is a 
unique material: the linear relaxation time for YIG is 100ns+. 
This allows dephasing of the wave packet in our experiment 
well before amplitudes of its constituents drop to the thermal 
level (see Fig. 3, where the typical time scale is hundreds of 
nanoseconds, which is given by the rate of dephasing.) For 
comparison, the best metallic magnetic material - Permalloy - 
has magnetic relaxation time < 5 ns. This material is opera- 
tional in the same frequency range as YIG, thus the dephasing 
rate would be the same, be this experiment conducted with 
Permalloy. Obviously, for the reason of the much stronger de- 
cay rate we would not be able to register any signal recovered 
by the second pump pulse, as on the time scale of hundred 
nanoseconds no coherent signal would survive in a Permalloy 
film to the time of its arrival. However, the new materials, like 
Heusler alloys and other half-metals, theoretically may have 
magnetic losses of the same order as YIG 28-30 . If quality of 
these materials is improved in the nearest future, it will be 
possible to use the two-pulse approach for these materials as 
well. 

V Conclusion 

In this work we have investigated the relaxation of a free 
evolving gas of previously parametrically pumped magnons. 
The experimental results show a clear deviation from the stan- 
dard exponential spin- wave decay model. In particular, the 
inherent magnon damping is found to depend upon the pres- 
ence of the parametric pumping field. The results are in agree- 
ment with the model which accounts for variation of phase 
coherence for parametrically injected magnon groups during 
the pump-free pause. 

Acknowledgment 

Financial support by the DFG within the SFB/TRR49 and 
by the Australian Research Council is gratefully acknowl- 
edged. The authors thank Alexy Karenowska for strong sup- 
port in the preparation of this paper. 



1 H. Suhl, IEEE Trans. Mag. 34, 1834 (1998). 

2 T. Gilbert, IEEE Trans. Mag. 40, 3443 (2004). 

3 A. Azevedo, A. B. Oliveira, F. M. de Aguiar, and S. M. Rezende, 
Phys. Rev. B 62, 5331 (2000). 

4 S. Rezende, A. Azevedo, M. Lucena, and F. de Aguiar, Phys. Rev. 
B 63, 214418 (2001). 

5 A. Y. Dobin and R. Victora, Phys. Rev. Lett. 90, 167203 (2003). 

6 K. Gilmore and M. D. Stiles, Phys. Rev. B 81, 174414 (2010). 

7 S. O. Demokritov, V. E. Demidov, O. Dzyapko, G. A. Melkov, 
A. A. Serga, B. Hillebrands, and A. N. Slavin, Nature 443, 430 
(2006). 

8 A. Khitun, M. Bao, and K. L. Wang, J. Phys. D: Appl. Phys. 43, 
264005(2007). 

9 O. Dzyapko, V. E. Demidov, S. O. Demokritov, G. A. Melkov, and 
A. N. Slavin, New Jour. Phys. 9, 64 (2007). 

10 H. Suhl, J. Phys. Chem. Solids 1, 209 (1957). 

11 R. M. White and M. Sparks, Phys. Rev. 130, 632 (1963). 

12 A. A. Serga, A. V. Chumak, A. Andre, G A. Melkov, A. N. Slavin, 



S. O. Demokritov, and B. Hillebrands, Phys. Rev. Lett. 99, 227202 
(2007). 

13 V. Cherepanov, I. Kolokolov, and V. L'vov, Phys. Rep. 229, 81 
(1993). 

14 Y. Kajiwara, K. Harii, S. Takahashi, J. Ohe, K. Uchida, 
M. Mizuguchi, H. Umezawa, H. Kawai, K. Ando, K. Takanashi, 
S. Maekawa, and E. Saitoh, Nature, 464, 262 (2010). 

15 K. Uchida, J. Xiao, H. Adachi, J. Ohe, S. Takahashi, J. Ieda, 
T. Ota, Y. Kajiwara, H. Umezawa, H. Kawai, G E. W Bauer, 
S. Maekawa, and E. Saitoh, Nature Mater. 9, 894 (2010). 

16 S. A. Manuilov, R. Fors, S. I. Khartsev, and A. M. Grishin, J. 
Appl. Phys. 105, 033917 (2009). 

17 O. V. Kolokoltsev, C. L. Ordonez-Romero, N. Qureshi, O. Cortes- 
Perez, G Lopez-Maldonado, and M. Avendano-Alejo, Electron. 
Lett. 46, 1387 (2010). 

18 R. W Damon and J. R. Eshbach, J. Appl. Phys. 31, 104 (1960). 

19 B. A. Kalinikos, N. G Kovshikov, M. P. Kostylev, and H. Benner, 
JETP Letters 64, 171 (1996). 



20 B. A. Kalinikos and M. P. Kostylev, IEEE Trans. Mag. 33, 3445 
(1997). 

21 S. Schafer, A. V. Chumak, A. A. Serga, G. A. Melkov, and 
B. Hillebrands, Appl. Phys. Lett. 92, 162514 (2008). 

22 A. V. Chumak, A. A. Serga, B. Hillebrands, G. A. Melkov, 
V. Tiberkevich, and A. N. Slavin, Phys. Rev. B 79, 014405 (2009). 

23 V. S. L'vov, Wave Turbulence Under Parametric Excitation 
(Springer, Berlin, 1994). 

24 V. E. Zakharov, V. S. L'vov, and S. S. Starobinets, Sov. Phys. Us- 
pekhi 17, 896 (1975). 

25 G. A. Melkov, Y. V. Kobljanskyj, A. A. Serga, V. S. Tiberkevich, 
and A. N. Slavin, Phys. Rev. Lett. 86, 4918 (2001). 



26 G. F. Herrmann, R. M. Hill, and D. E. Kaplan, Phys. Rev. B 2, 
2587 (1970). 

27 M. M. Scott, C. E. Patton, M. P. Kostylev, and B. A. Kalinikos, J. 
Appl. Phys. 95, 6294 (2004). 

28 T. Kubota, S. Tsunegi, M. Oogane, S. Mizukami, T. Miyazaki, 
H. Naganuma, and Y. Ando, Appl. Phys. Lett. 94, 122504 (2009). 

29 C. Liu, C. K. A. Mewes, M. Chshiev, T. Mewes, and W. H. Butler, 
Appl. Phys. Lett. 95, 022509 (2009). 

30 S. Trudel, O. Gaier, J. Hamrle, and B. Hillebrands, J. Phys. D: 
Appl. Phys. 43, 193001 (2010). 



